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ABSTRACT 


The intent of this paper is to describe a method of coloring a 
map and to present an algorithm for the solution of this problen. 
A camputer program was developed to provide solutions to the problem 
of coloring a map which consists of a finite number of areas. This 


algorithm may also be applied to problems other than map-coloring. 
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SYMBOL 


TABLE OF SYMBOLS 


MEANING 


Number of equations used to describe the border 
constraints for a map. 


Number of areas to be colored 

Variable representing the color of area 1 

The MXN matrix defining regions having common borders 
The 5) column of Ayn 

The i row of Avan 

ee la, ee ee 


The number of zero elements of each colum a. that 
are not associated with the constraint equations 


The number of zero elements of each column a. 





I. INTRODUCTION 


When coloring a geographical map, it is custamary to use different 
colors for adjacent regions - where adjacent means a canmon finite 
boundary. It has been conjectured, but not proved, that "For any sub- 
division of the plane into non-overlapping regions, it is always pos- 
sible to mark the regions with one of the numbers 0, 1, 2, 3, in such 
a way that no two adjacent regions receive the same nurber "+ 

The statement that a map can be colored with four colors is 
accredited to Moebius, who first proposed it in 1840. This conjecture 
has never been proven; however, a satisfactory counter-example has not 
been demonstrated either. A proof for five-coloring a map may be found 
in [Courant and Robbins 1941] and [Oystein 1967]. 

This thesis will use the four-coloring conjecture without proof 
and present an algorithm that utilizes an iterative routine to determine 
what color each region should be. A similar iterative technique for 
general integer programming problems is demonstrated in [Greenberg April 
1969] and [Greenberg May 1969], which present the proofs. We are 
interested in showing whether the iterative technique can be used for 
large-scale problems such as the map-coloring problem. 

The problem formulation suggested by Gomory [Dantzig 1960] is used. 
Let the map regions be denoted by i = 1,2,...,N and let Xs be an integer 


valued variable such that 


Oe Ss (1) 


Bt ae, Richard and Robbins, Herbert; What is Mathematics?, p. 247, 
Oxford University Press, 1941. 


where the four values x. = 0, 1, 2, or 3 correspond to four different 
colors. Since any two adjacent areas must have different colors, the 


constraint may be written in an either/or form 


either x. —- xX. 3°1 = - xX. SL (2) 
i | 5 i 


if i and j have a common border. Equation (2) may be rewritten as: 


Xx. - X. > 1 - 46.. (oe = 0) 
a J 1) 1) 


(3) 
Xx. - X. > 3 + 46... 
J ~ 1) 


Instead of using (2) and (3) we can simplify the problem by using 
ie x 0 (4) 


where 1 and j have a cammon border. 
This formulation can be written in matrix form 

Ax # 0 (5) 
where the elements of Ava aeu are either 0, -l, or 1. We can now 
enumerate all values of the left side of (5). The column with the 
-least number of zeroes is now selected from Aya a, and the value 
of the corresponding x, is set equal to one. The colum a. is then 
added column by column to the matrix Ay to form a new matrix. The 
process continues in the same way by searching the new matrix for the 
column with the least number of zeros. In each iteration the correspond- 
ing = is increased by one. The column values represent the value of 
the left side of (5). A solution occurs when a colum is achieved that 


has no zeros. 


“The Speed with which the algorithm reaches a solution is increased 
by modifying the method of choosing a.. This modification is explained 
in Section II. A. 4. : 
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The algorithm presented is based upon the four-color conjecture, 
but is not restricted to the use of only four colors, if more colors 


are desired. 


II. THE COLORING PROCESS 


A. THE ALGORITHM 

1. Let x, = 0 for 1 = 1,...,N. 

2. For every 2 regions of the map with a camon border it is 
required that a 70, 1d = Lo dll ¥~ 5. 

3. Form the M x N matrix (Any defined in step 2, where M is 
the number of equations required to describe the border constraints 
and N is the number of areas to be colored. The elements of Aven Ay 5! 
are either -]1, 0, or 1. Note: 

a. For each row, a., 1=1,.-.,M, there are only two non zero 
elements, -1, and 1, corresponding to the constraints in step 2. 

b. For colum ay, i =1,..06N, has at least one monmzene 
element and fewer than M non zero elements. 

4. Let ae j =1,...,N, be the number of zero elements of each 
a.. Let 8., j = 1,...,N, be the number of zero elements of each a, 


zy i) 


which appear in those positions where ay5 = QO. Note: B < Cig 


g = 1,...,N. Find the coli of Ayn with the smallest BS Let ae* 
be the smallest j such that Baas = min iB. }. het 9°" ==. 
5. Let x'., = x., + 1. Check all the constraints in which xT as 


=| 
appears to insure that the constraint equations (4) are satisfied. If 
the constraint equations (4) are not satisfied, continue increasing 


Xi in increments of one until either (4) is satisfied or ‘oe = 3. 


Tf ee = 3 and (4) is not satisfied, set ae = 0 and let j* be the 


dele 


smallest j > j** such that 8., = min {8.3 and repeat step 5. If no 


ae 

such j* exists, then find the — a and proceed as with the 

BS. If no ae exists such that Xo and x. satisfy (4) for all i # j*, 
go to step 9. 

6. Determine the elements of column a.,, by substituting the presently 
defined values of Xe, i =1,...,N, into the constraint equations (4). 

Add the column Bsr j =1,...,N, to each colum of A... Let A’ = 
ee Lasers es ss yl 

Tee ee Oe j = 1,...,N, be the number of zero elements of each a's. 
Let Bs, j = 1,...,N, be the number of zero elements of each a’ which 
appear in those positions where mae = 0. Let j** be the smallest j such 
that be = (8.3 o Sta) sae. 

8. Continue steps 5 through 8 until = = 0, in which case a solution 
has been reached. Note: If x. = 3 at any point during steps 5 through 
8, disregard the corresponding column a. in any future iteration, since 
the corresponding x. is at its maximum possible value. 


MxN 
matrix. For each zero element of ass there are 2 regions, g and h, with 


9. Let the colum as x be the one used to generate the latest A! 


a common border where e secre 0, 1.e., a Violation of (4). Let F be 
the set of all such regions g and h. Denote by the subscript f the 
element of F. 


10. For each f,increment the x,'s such that Xe = Xp + 1. Check all 


i 


of the constraints in which xX, appears to see if the constraint equations 
(4) are satisfied. If for any f, an element of F, the constraint equations 
(4) are satisfied set f = j* and go to step 6. If (4) is not satisfied 


for any i, not an element of F, where regions i and f have a cammon 


M2. 


Initial 




















Symbols: 
PS (1) in direction of 
> decreasing we 
--— - > back tracking 
am 80785 (9) ] 
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J(f) 
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Possible Path to a solution to the map-coloring problem 


Figure l 


Migs 


border, note the iteration that defined the present value of X.. 
Repeat step 10 until Xe = 3. Reset Xe = Q. Repeat step 10 for all 
He 
ll. Order the iterations defined in step 10 in decreasing 
iteration number. Let these iterations be called possible iterations. 
12. Reconstruct the most recent possible iteration and examine 
B = Bas Se ne j = j* + 1,...,N, depending on which was used 
to determine Xa in the iteration being examined. Choose a new j* 
in the same manner as the j* was chosen in step 5. When a j* is found 
such that Xa and x. satisfy (4) for all i # j*, go to step 6. If no 
such j* exists satisfying (4), then the iteration being examined is not 
possible. Repeat step 12. 
This completes the algorithm. If no possible solution exists, the 
algorithm will continue "backing up" the generated matrix chain until 
the original matrix is reached. Figure 1 shows a possible path to the 


solution of the map-coloring problem. 


B. EXAMPLE 





Fig. 2. Example Problem 


>the speed with which the algorithm reaches a solution is increased 
by starting the backtracking procedure when no variable associated with 
the minimum B. or as can be set at any integer value between 1 and 3. 
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where 


Xx) — X5 x 0 
x) X x~ 0 
x) - Xy # 0 


set x. = 0 forall 


we write 


column used * 


then xy = I 


the column to add (a. 


5%? 


ie 
7 antans 
2 | = vA 


where as = [111000] 


then adding Cone 5 


Xo — Xy x 0 
X — Xy x 0 
X. ~ X, x 0 
X, XX, xy 
=I. 0 0 
Opes 0 
0 0 -)j 
1 Si 0 
i elt 
0 el 
3 3 3 
3 3 5 


1s defined by: 


i Xo ~ x, = 0 
l Xo ~ xy = Q 
1 x, ~ xy = 0 
+ Anan we write 


nS) 


column used * 


then Xx. = 2 


the column to add (a5) is defined by: 


xX ~ X> = = Xo» ~ 
xX) — X3 = 1 x 
Xx) — X4 = J, x3 - 


where Bye = [-111220] 


then adding lay gr ++ rsx! + Ayn We ee 


1 -l 
uae 1 1 
J 
8. 1 0 
column used “ 
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then X., = 3 


the colum to add (Aj) is defined by: 


x) ~ Xo = -] Xo - Xo atl 
x) = x, = -2 x. ~ xy = 7, 
x) - xX, = 1 X3 oot 3 


where :... = [ee Se 


then adding [ajare ++ 1854] +A, we write 


eee eo oe 
0 =) il 
= =) eee? 
2 1 1 
= le SB 
z 3 2 
3 3 4 
ro) 1 ue 0 


there is no need to calculate column a', because there are no zeros 





4 
in colum a’, and a solution has been reached. The solution is: 
Area number Color number 
Il al 
Z, 2 
3 3 
4 0 


Wy 


IfifL. CONCLUSION 


The algorithm described will solve the map-coloring problem. 
Although only four colors were used in the example in Section II-B, 
the algorithm can be adapted to solve a problem for any feasible 
number of colors. Several representative map-coloring problems were 
solved using the enclosed computer program. Appendix A shows a 
solution to four-coloring the continental United States. In the 
representative problems solved, both four and five colors were used. 

If the solution is not feasible, e.g., the use of four colors 
when five colors are required, as in the split-state problem, the 
algorithm may consume several hours of computer time before showing 
that no solution exists. The computer time is consumed in back- 
tracking the solution path and in investigating every possible 


solution path for both the a's and the 8's. 


18 


APPENDIX A 


SAMPLE PROBLEM SOLUTION 


The solution for four-coloring the continental United States of 
America, including boundaries and major water bodies, is as follows: 
STATE COLOR 


Maine 

New Hampshire 
Vermont 
Massachusetts 
Rhode Island 
Connecticut 
New York 

New Jersey 
Pennsylvania 
Delaware 
Maryland 
Virginia 
North Carolina 
South Carolina 
Georgia 
Florida 
Alabama 
Tennessee 
West Virginia 
Ohio 

Indiana 
Kentucky 
Illinois 
Wisconsin 
Michigan 
Minnesota 
Iowa 
Missouri 
Arkansas 
Mississippl 
Louisiana 
Texas 

Ok lLahoma 
Kansas 
Nebraska 
South Dakota 
North Dakota 
Montana 


Wyoming 


NO OWFOWONONWNHOWOWNDOWNOKRFEFWONOWNOWNHOWOWN OC WO 


ANS, 


Colorado 
New Mexico 
Arizona 
Utah 

Idaho 
Washington 
Oregon 
Nevada 
California 


WONM OF WF WF 


BORDER AREAS 


Canada 

Mexico 

Water - including: 
Atlantic Ocean 
Pacific Ocean 
Great Lakes 
Great Salt Lake 


"oS 


The above problem was solved in 97.2 seconds on the IBM 360/67 camputer 


and used 81 iterations. 
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Yes 


APPENDIX B 


PROGRAM FLOW CHART 


Assign initial 
values to vari- 


ables 





Read in parameter 
and border 
constraints 










Calculate the # 
of cammon borders 
for each area 







Calculate aj, 63 
j =1,...,N and 


Avan a Oe 






zero the counter 
fer dak ee, 
4 =1,.2.,N 





No 


Count iteration 
number 


| Calculate BX 


Store a.,/;X.4/)* 
check x., 


Mark the column 
with the smallest 


B. (Denote as j*) 
Is 
there a 
possible 
solution? 





No 
STOP 
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peeme da x. m,)* 


Return 


LJ 
LI 


Yes 
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there a 
possible 
solution? 







Yes 


Determine which 
variables still 
= 0 and have com- 
mon borders 







For above vari- 
ables set 


Xe = Xe + ] 








Calculate which 
lteration defined 
value of Xe 









variables been 
examined for all 
possible colors? 







Yes 
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No 


No 


Return 


Define a new_ 
column ae = 2. 


Order the iter- 
ations in de- 
creasing iteratic 


ii 






Return to most 
recent possible 
iteration 








Yes 





Assign to each 
variable the value 
used prior to mos 
recent possible 

iteration 








reconstruct chosen 
iteration 


/ 


Use next column 
with same number 
of zeros 











If no colum exists 
with same number 
of zeros, return 
to next possible 
iteration 
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No _: 


No solution exists 
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